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ABSTRACT 

Non-linear nature of Einstein equation introduces genuine relativistic higher order corrections to the usual 
Newtonian fluid equations describing the evolution of cosmological perturbations. We study the effect of such 
novel non-linearities on the next-to-leading order matter and velocity power spectra for the case of pressureless, 
irrotational fluid in a flat Friedmann background. We find that pure general relativistic corrections are negligibly 
small over all scales. Our result guarantees that, in the current paradigm of standard cosmology, one can safely 
use Newtonian cosmology even in non-linear regimes. 
Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

Large scale structure (LSS) of the universe is a powerful 
probe to study the nature of cosmological d ensity perturba - 
tions and to extract cosmological parameters (Peebles 1980). 
Combined with the anisotropy of the cosmic microwave back- 
ground (CMB), most of the cosmological parameters are cur- 
ren tly constrained within a few percent accuracy or even bet- 
ter dKomatsu et al. 2010). To continue our success in cosmol- 
ogy with CMB and LSS, it is crucial to predict the power spec- 
tra from theory accurately. While the temperature fluctuation 
in the CMB is as small as ST/T ~ 10" 5 dSmoot et alJI 19921) 
so that linear perturbation theory is able to provide necessary 
accuracy, we have larger degree of non-linearities in LSS. 
We must take into account non-linearities of LSS properly 
to predict the power spectrum accurate en ough for precision 
cosmology at a level similar to the CMB dJeong & Komatsul 
20Ml2009h . 

Most studies on LSS, however, have been based on New- 
tonian gravity, especially those including non-linear perturba- 
tion s dVishniac|[T983l:lGoroff et all 1986: M akino et al.llT992i: 
lFry||1994[ iBernarde au et alj|2002l) . This approach has to be 
justified a posteriori by comparing th e result against fully g en- 
eral relativistic one. For example, in lNoh & Hwangl (12004). it 
is shown that the Newtonian hydrodynamic equations up to 
second order coincide exactly with the relativistic ones in the 
zero pressure case, after appropriately identifying hydrody- 
namical variables with gauge-invariant combinations of rela- 
tivistic perturbation variables. Thus, compared to the Newto- 
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nian hydrodynamic equations which are closed at second or- 
der, any higher order contri butions are originated from purely 
general relativistic effects (Hwang & Noh 2005b). A consis- 
tent expansion of density fluctuation tells us that the leading 
non-linear contributi ons to the power spec t rum include thir d 
order perturbations dNoh & Hwangl 120081: iNoh et alj|2009h . 
Thus non-linear density power spectrum naturally include 
pure general relativistic effects, which may have important 
implications as N-body simulations are becomin g larger and 
larger to reach the horizon scale (Kim et al. 20Q9|)- 

In this note, we examine the general relativistic effects on 
the power spectra of matter density fluctuations and peculiar 
velocity by including leading non-vanishing non-linear con- 
tributions. Our aim is to answer the question whether pure 
general relativistic effects can give rise to any cosmologically 
observable consequences. To our surprise, we find that the 
Newtonian terms in these power spectra are absolutely domi- 
nating over all relevant cosmological scales, even outside the 
horizon. Although the result sounds simple and pleasant, this 
is still a non-trivial result, because in the context of cosmol- 
ogy Newtonian gravity is incomplete: there is no concept of 
horizon, the propagation speed of an action at one point is 
infinite, and so on. 

This note is outlined as follows. In Section |2j we present 
the formalism to set up the equations to solve, and give the 
solutions up to third order. In Section[3j we compute the mat- 
ter and velocity power spectra including next-to-leading non- 
linear corrections which include genuine general relativistic 
effects. In Section|4]we conclude. 

2. EQUATIONS AND SOLUTIONS 

We consider Einstein-de Sitter universe, i.e. a flat universe 
dominated by pressureless, irrotational matter, and consider 
only the scalar perturbations. We work in the temporal co- 
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moving gauge where T°i = to all perturbation orders, with 
i being a spatial index. As this temporal gauge c ondition, 
together with our unique spatial gauge condition (Bardeen 
1988) gij = a 2 (l + 2(p)8ij, fixes the gauge degrees of free- 
dom completely, all the resulting perturbation variables can be 
equivalently regarded as fully gauge invariant, both spatially 
and t emporally. This state ment is valid in all perturbation or- 
ders (lNoh & Hwandl2004l) . 

The Arnowitt-Deser-Misner formulation ( Arnowitt et all 
2008) is convenient in our case (Bardeen 1980). The comov- 
ing gauge condition imply the momentum density vanishes, 
i.e. Ji = NT°i = 0, with N being the lapse function. The pres- 
sureless condition implies Sy = 7y = 0. Therefore the mo- 
mentum conservation equation gives Nj = 0, thus the lapse 
function N is uniform. The energy and momentum conserva- 
tion equation s and the trace part of the propagation equation 
then become dBarde en 1980) 

Efl-N'Ej = NKE , (1) 
^-!*|i = 0, (2) 

K fi - N'Kj = n(^K 2 + K l jK J i + 4nGE - A^j , (3) 

where Nj is the shift vector, K = K', is the trace of the ex- 
trinsic curvature tensor K,j, A is the cosmological constant, 
£=Af 2 r 00 isthe energy density, an overbar denotes the trace- 
less part, and a vertical bar denotes a covariant derivative with 
respect to Jh. These are the complete equations we need in 
our non-linear perturbations, and are valid in fully non-linear 
situation. 

We introduce the density and the velocity fluctuations as 
E = p(t) + 5p(t,x) and K= 3H-6(t,x) with 6(t,x) = a"'V • 
u(t,x), with a being the cosmic scale factor. We can iden- 
tify Sp(t,x) and u(t,x) as the Newtonian density and veloc- 
ity perturbation variables respectively, because the relativistic 
equations coincide exactly with the corresponding Newtonian 
hydrodynamic equations up to second order. 

From the above equations we can derive the hydrody- 
namic equations of density fluctuation 8(t,x) = 5p(t,x)/ p(t) 
and ve locity gradient 9{t,x) to the third order (Hwan g & Nohl 
2005b). The relativistic continuity and Euler equations are 
found to be 

^ + -V-lt = --V-(&») + - \2ipu-V (A- l X 2 )] -(V<5), 
at a a a 



(4) 



- V ■ ( ^-+Hu) + 4nGpS = ~ V • [(« • V)«] 
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(« • V)m- — (V • u)u 
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3a 



z ip(u- V)(V-w) 



- 1a [(« ■ V)A~ l X 2 ] + 1(b • V)I 2 + ^I 2 (V • u) , (5) 

where (p and X 2 are the linear and the second order quantity 
respectively, and are defined as 



In relativistic perturbation theory, the dimensionless quantity 
ip is proportional to the spatial curvature perturbation in the 
comoving gauge. All the perturbation variables 5, u and tp 
can be regarded as equivalently gauge-invariant combinations 
to non-linear order. Proper choice of variables and gauge 
conditions are important to have these equations. Note that 
the relativistic continuity and Euler equations coincide with 
those from Newtonian fluid approximation up to the secon d 
order in perturbations (Peebles| [l980t lNoh & Hwangl l2004). 
Therefore, the perturbative solutions are also the same up to 
the second order, and pure general relativistic effects appear 
from third order. We emphasize that the above equations are 
valid in the presence of the cosmological constant in the back- 
ground world model. 

An examination of the third order terms in Equations (|4|i 
and (O shows that the pure third order terms are simple con- 
volutions of the linear order <p with the second-order com- 
binations of fluid variables S and u. Note that to the linear 
order, <p is a well-known conserved quantity whose ampli- 
tude of the growing mode solution is conserved on the super- 
horizon scales, independent of the changing e quation of state 
or even changing underlying gravity theories dHwang & Nohl 
I2005al) . In a flat background without cosmological constant 
the amplitude of <p near horizon scale is directly related to the 
amplitude of relative temperature fluctuations of the CMB as 
8T/T = ip/5. 

The linear solutions of Equations (|4|i and (|5]l are easily 
found to be 



6 1 (k,t)=D(t)Si(k,t ), 
6 l (k,t)=-aHD(t)5 l (k,t ), 



(8) 
(9) 



where D(t) is the linear growth factor so that 8\(k,to) is the 
present linear density fluctuation. With these linear solutions, 
we can perturbatively expand the density contrast 5{k,t) = 

61 + 62 + 83-1 , where 8„ is a n-th order quantity in linear 

density contrast 6i(k,to), and similarly for 6(k,t). With this 
expansion, we can find the full non-linear solutions of Equa- 
tions (H]i and (|5j by using momentum dependent symmetric 
kernels as 



n- 1 \ i=l / 

*F^(q u ---q n )6 l (q l )---5 l {q n ), (10) 
8(k, t) =-aH±D» (t) J (k-±q) 



xG^(q 1 ,-q n )6 1 (q 1 )-6 1 (q n ). 



(11) 



^tp = i ( -\ixGpb+ — V ■ u ) . 



(6) 



X 2 = 2ipV -u-(u- V)ip+-A~ l V -[uAip+u-V(Vtp)] . (7) 



Then, Equations (0|i and (0 become simple differential equa- 
tions of and G^K Especially, the general relativistic terms 
which explicitly include kjj = aH, the comoving wavenumber 
corresponding to the comoving horizon, are reduced to the al- 
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where we have introduced # 12 ...„ = Y^i=i9r The second and 
third order Newtonian kernels can be fo und in e.g . Equations 
[(2.32), (2.33)] and [(2.34), (2.35)] in Ueong I d2010h . respec- 
tively. 

3. MATTER AND VELOCITY POWER SPECTRA 

From Equation ( [Tol l, we can find the non-linear power spec- 
trum, which is defined as 

(5(k u t)5(k 2 ,t)) = (27rfS 0, (k 1 +k 2 )P(k u t). (14) 

If we assume perfect Gaussianity of 8\, which is a very 
good approximation consistent with current observations, any 
higher order correlation function beyond the linear power 
spectrum Pi\(k) disappears and Pn(fc) is all that we need to 
specify the statistics of density fluctuation 5: as we will see 
shortly, all the non-linear corrections to the power spectrum 
can be written in terms of P u . Then, from Equation (TT4l we 
can write, beyond the linear density power spectrum Pn, 

P = Pn+Pi2 + Pn + --- , (15) 

with (Si(ki)Sj(k 2 )) = (2TT) 3 S (3) (k l +k2)s i jP i j(k l ). Here, s u is 
a symmetric factor which is 1 for i = j and 1 /2 otherwise. 
The leading non-linear correction P\2 includes bispectrum and 
thus disappears according to our assumption of Gaussianity 
of 5\. P22 + P13 denotes the next-to-leading order non-linear 
correction to the power spectrum. As mentioned above, Pi 3 
includes general relativistic terms. 

The density power spectrum up to next-to-leading order 
non-linear corrections is 



P(k,t) = P n (k,t) + 
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=Pll +P22 + Pl3,Newton + Pl3,Einstein ■ 



1 -r 
(16) 



where r and x are the magnitude of dummy integration mo- 
mentum q and the cosine between q and k, respectively, intro- 
duced as q = rk (0 < r < 00) and k q = k 2 rx (-1 < x < 1). 
We have divided P13 into the Newtonian part P^Newton and 
the general relativistic contribution Pn.Einstein- Compared with 
Pi3,Newton me general relativistic contribution P[ 3 Einstein is mul- 
tiplied by a factor {kn/k) 2 , where ka/k is the ratio between a 
scale of interest and the horizon scale, and is thus highly sup- 
pressed far inside the horizon. 

In Figure Q] we present the total power spectrum of Equa- 
tion dT6l > along with its components Pn, P22, Pi3,Newton 
and Pi3.Einstein when our Universe is dominated by mat- 
ter, at z = 6. The l i near p ower spectrum is calculated by 
CAMB dLewis et all 120001) code with the maximum like- 
lihood cosmological parameters given in the Table 1 of 
iKomatsuet ail (2009) ("WMAP+BAO+SN"). FigureQ]shows 
that the general relativistic contribution Pi3,Einstein is smaller 
than the linear power spectrum P\ 1 on all cosmological scales. 
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Figure 1. Non-linear matter power spectr um ( solid blue line) and the con- 
tribution from each component of Equation (16) at z = 6. The black, red, and 
orange lines show the contributions from the Newtonian perturbation the- 
ory: Pu, P22, and Pi3.Newtom respectively. The green line shows the general 
relativistic effect, f 13. Einstein W- Note that we take the absolute value for neg- 
ative terms, and show with dashed lines: P22 an d Pn are positive and neg- 
ative,respectively, in all scales. Vertical dotted line shows the wavenumber 
correspond to the comoving horizon kn at z = 6. Over all scales, the general 
relativistic term Pi3.Einstein (green) is negligibly small compare to the linear 
power spectrum P\ \ (black). 

Let us examine Pn.Einstein more closely. For notational sim- 
plicity, we shall abbreviate the integration in Pi3,Einstein as 
J drP\i(kr,t)f(r). Then, scale dependence of Pi3,Einstein can 
be understood as follows. First, setting kr = q, we find that 
Pn.Einstein ~ P\\(k)f(q/k). On small scales (k > 0.01/i/Mpc), 
q/k is also small, and by using Taylor expansion of f(r) = 
-(656/15)r 2 +C(r 4 ) we find Pi 3 , E instein ~ k~ 2 P n . On the other 
hands, in large scale limit (k <C 0.01/i/Mpc) where q/k takes 
larger value, f(r) = -752/3 + (D(r~ 2 ) and P^Einstein has a scale 
dependence Pn.Einstein ~ Pn- Numerical calculation reveals 
that Pi3,Einstein is smaller than Pn by a factor 10~ 5 on large 
scales. Our result shows that the leading order non-linear 
power spectrum is finite in both infrared and ultraviolet re- 
gionfl 

2 The previous result reporting infrared divergence in Pi3.Einstein 

INoh et alj 
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We can proceed in the same way to compute the power 
spectrum of the peculiar velocity. As Equation ( TBI ), we can 
define 



(9(k u t)6(k 2 ,t)) = (2ir) 3 S (3 \k l+ k 2 )P ee (k u t), 
and we can find 

1 k 3 



(17) 



k- H l P gg (k,t) = P n (k,t) + 
• i 



drP u (kr,t) 



dxPu (k^l + r 2 - 2rx, t 



98 (2tt) 2 jo 

{r-lx+brx 2 ) 1 
(\ + r 2 -2rx) 2 
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+ ^7T-a Pn ^ k '^ / drp n(kr,t) 
84 (27r) z Jo 
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l+r 
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(18) 



Figure [2] shows the non-linear velocity power spectrum of 
Equation (fT8l ) for exactly the same cosmology as Figure [1] 
As in the case of the total matter power spectrum, the non- 
linear general relativistic correction is negligibly small for 
all scales. It is because the third order kernel for velocity 
(^Einstein behaves in the same way as that for the matter den- 
sity /^Einstein in both large (r — > 0) and small (r — > oo) scale 
limit: km^ogCr) = -(368/15)r 2 + C(r 4 ) and linWoo g(r) = 
-496/3 + Oir~ 2 ) when denoting the last integration in Equa- 
tion ylO as / drP n (kr)g(r). 
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Figure 2. Same as Figure[T] but for the velocity power spectrum Pgg(k). 



To conclude, in this note we have examined the general rel- 
ativistic non-linear contributions to the density and velocity 
power spectra. We have found that, with pleasant surprise, 
the pure general relativistic effects are completely negligible 
on all cosmologically relevant scales, even outside the hori- 
zon. It is interesting to see that the linear power spectrum is 
totally dominating even outside the horizon. Our conclusion 
has the following important implication. As the general rel- 
ativistic effect is very small, Newtonian theory can be safely 
applied to the non-linear evolution of cosmic structure on all 
cosmologically relevant scales. In the literature it has been 
common to use Newtonian gravity to study the non-linear 
clustering properties of large scale structure without justifying 
that approach. The result we present in this note provides a 
confirmation of using Newtonian gravity to handle non-linear 
clustering in cosmology. 
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4. CONCLUSION 



2009) turns out to be due to an incorrect calculation of the power spectrum: 
the third order general relativistic kernels F3 Einstein nas not been fully sym- 
metrized, thus causing logarithmic infrared divergence in Pi3.Einstein- 



